Polydispersity-linked Memory Effects in a Magnetic Nanoparticle System 
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We have performed a series of measurements on the low temperature behavior of a magnetic nano- 
particle system. Our results show striking memory effects in the dc magnetization. Dipolar in- 
teractions among the nano-particles suppress the memory effect. We explain this phenomenon by 
the superposition of different super paramagnetic relaxation times of single domain magnetic nano- 
particles. Moreover, we observe a crossover in the temperature dependence of coercivity. We show 
that a dilute dispersion of particles with a flat size distribution yields the best memory. 

PACS numbers: 75.75.+a, 75.50.Lk, 75.50.Tt, 75.20.-g 
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Single domain magnetic nano-particlesQ, 0, El 0, IE IE 
„ possess relaxation times that depend exponentially 
on the volume H El El E2- Thus polydispersity leads 
to a distribution of relaxation times|ia.ll4|. those larger 
than the measurement time yielding 'frozen' behavior, 
and those shorter giving rise to 'magnetic viscosity' [ill 
If3 |. A given sample then displays strong memory effects, 
which are reported here. Our results are based on the 
measurements of temperature-dependent magnetization 
during cooling and heating cycles. These memory effects 
may have important device applications^ in the future. 

In this Letter we report results on magnetization and 
coercivity measurements in systems of NiFeiO^ parti- 
cles (mean size~ 3nm.) embedded in a SiOi matrix. 
Both measurements show strong history dependent ef- 
fects depending on the separation between the particles 
and hence their mutual interaction. We prepare the 
samples by using a sol-gelE3 technique. The ratio of 
NiFeiO^ to S1O2 is 1:1 and 3:7 by weight, for sample A 
and sample B respectively. The ferrite to marix ratio is 
altered in order to tune the interparticle interaction. The 
phase of the samples is identified by X-ray diffractionfl^ 
using a Philips PW 1710 diffractometer with Cu K Q ra- 
diation (A = 1.54 A°). The average particle size is es- 
timated by X-ray diffraction as well as by a JEM-200- 
CX transmission electron micrograph(TEM). DC mag- 
netization measurements are performed on a Quantum 
Design superconducting interference device magnetome- 
ter (MPMS) from 300k to 4k. 

The X-ray diffraction spectrum confirms that our sam- 
ple is indeed in a single-phase of ferrite NiFe2 04 with 
no residual a-Fe203. The average particle size of the 
nanophase NiFe2 04, estimated from the broadening of 
X-ray diffraction peak is ~ 3nm for both the samples. 
The TEM micrograph of the samples suggests that the 
particles are spherical in shape and the particle size fol- 
lows a log-normal 14J distribution. The interparticle sep- 
aration measured from TEM micrographs are ~ 5nm for 
sample A and ~ 15nm for sample B. 
The magnetization measurements are carried out in ac- 



cordance with the following cooling and heating proto- 
col. At T=300K (T = Too), a small magnetic field 
(H = 50 Oe) is applied and the magnetization(M) mea- 
sured. Keeping the field on, the temperature (T) is low- 
ered continuously at a steady rate to T n and M is si- 
multaneously measured upto temperature T n . Thus one 
obtains M versus T in the cooling regime (T n <T< Too). 
At T„ the field is switched off and the drop of M is 
monitored for several (« 4) hours. Subsequently, the 
magnetic field is switched back on and M(T) versus T 
is mapped in the cooling regime(T n -i < T < T n ). At 
T„_i the field is switched off again and the process of 
measurement repeated, until the lowest temperature To 
is reached. Thus, one obtains field-cooled response and 
zero-field relaxation of the magnetization as a function 
of temperature. At the end of the cooling cycle, at T , 
the field is turned on and M(T) monitored as the system 
is heated from T through T n _ 2 ,T„_ 1 ,T„ and eventu- 
ally to Tqo, the magnetic field remaining on throughout. 
Our results are shown in Fig.l, for two distinct values of 
interparticle spacing 5nm. (sample A), and 15nm (sam- 
ple B). The heating path surprisingly shows wiggles in 
M(T) at all the T steps T n _2,T n _i,T„ where H was ear- 
lier switched off during cooling, apparently retaining a 
memory of the temperature steps at which the cooling 
was arrested. 

One tantalizing aspect of our results is that memory ef- 
fects are more prominent for sample B than for sample A, 
although in the latter the average inter-particle distance 
is smaller and hence the dipolar interaction larger. Re- 
cently Sun et. q/.|l3| have reported very similar history 
dependent effects in the magnetization measurements of 
a monolayer of sputtered permalloy(-ZVi 8 i.Fei9) clusters 
on a S1O2 substrate. These authors attribute the dis- 
parate cooling and heating histories to aging and con- 
comitant memory effects found in a spin glass phase E3- 
Spin glass transitions are known to occur due to dis- 
order and frustration^^ i n dilute magnetic alloys that 
are characterized by a complicated free energy landscape 
with deep valleys and barriers. Strongly nonequilibrium 
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FIG. 1: (color on-line) Experimental M(T) curves during cool- 
ing (□ red) in a small magnetic field H = 50 Oe and zero- 
field heating (O black) for the (a) interacting and (b) non- 
interacting samples showing prominent memory effects. A 
constant heating/cooling rate of 2° K was maintained except 
at 60°, 40° and 20° K where the cooling was arrested for 
4 h duration at each temperature during which time H was 
switched off. 



memory dependent behavior ensues as a result of the sys- 
tem getting trapped in a deep valley^] such that the re- 
laxation time(r) for deactivation becomes long compared 
to experimental time scales of measurement. 
Our interpretation of the results shown in Fig^ is very 
different from that of 01 • We demonstrate below 
that the observed phenomena (by us as well as by other 
groups) are not connected to complicated spin glass type 
interactions but can be simply attributed to a superposi- 
tion of relaxation times, arising from particle size distri- 
bution, as it were in noninteracting single-domain mag- 
netic particles. Experimentally it is known [l^ that nano 
particle sizes are usually distributed according to a log- 
normal distribution. However, we show below that the 
exact form of the distribution is irrelevant for explaining 
the memory effect. In fact, in order to keep the analysis 
simple and to obtain a clear understanding of the physics 
it is sufficient to take a simple size distribution consist- 
ing of two delta function peaks so that there are only 
two kinds of particles "large" (volume Vi) and "small" 
(V2). Correspondingly we have only two relaxation times 
t = n and r 2 in our model, if we remember 

nana 

that: 



t(V) oc exp 



(KV ± iiHV) 
Kb~T : 



(1) 



where K is the anisotropy energy, /1 is the magnetic mo- 
ment per unit volume, H is the applied magnetic field 
and Kb is the Boltzmann constant O ur interpre- 
tation of the observed results hinges on the premise 




FIG. 2: (color on-line) Simulated M(arbitrary unit) vs. T 
curves during cooling and heating for the (a) interacting and 

(b) non-interacting cases: The curve (c) shows the various 
contributions to the total magnetization (thick solid line)of 
interacting sample(A) coming from the fast (thin solid line) 
and slow (dashed line) particles. The theoretical curves (a) - 

(c) have been calculated using a double delta function distri- 
bution of particle sizes. Curve (d) shows a plot of the recov- 
ery parameter R (see text) as a function of the width(s) of a 
Gaussian particle size distribution. 



that the time T\ is much larger than the measurement 
time while r 2 is much smaller, at the lowest measured 
temperature (To). Both t\ and r 2 are expected to 
be smaller than the measurement time at the highest 
temperature Too. Therefore, in the intermediate tem- 
perature domain (To < T < Too), the small parti- 
cles equilibrate rapidly, thus showing superparamagnetic 
viscosity ^3 while the large particles are 'blocked'. This 
is observed in Fig^c) where we have plotted computer 
simulations [20j of M(T) separately for the two sets of 
interacting(separation distance of ~5nm) particles un- 
der the same cooling and heating regimens. Here we 
choose the temperature T* at which H is switched off 
such that the blocking temperatures^, ^| correspond- 
ing to the two different particle sizes flank T*.The simu- 
lations are based on standard rate theory calculation for 
the time dependent magnetization J2jj. When H is zero, 
both sets of particles relax to M=0. However, when H 
is turned on, particles 1 are blocked(M=0) while 2 show 
facile response. As T is increased again, M for particles 2 
decreases with T while M for particles 1 initially increases 
before dropping off. The resultant graph is a superposi- 
tion (sec FigEl (c)) of a monotonically decreasing curve 
and a hump, thus producing a wiggle. This effect is seen 
only when the temperature of arrest is in-between the two 
respective blocking temperatures, in conformity with the 
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findings of [ldj . 

We have performed measurements on the same sys- 
tem but now with increased interparticle separation 
(~15nm.)(see Figmb)),the simulation results of which 
are shown in Figl^b). 

The resultant interaction effect due to dipole-dipole cou- 
pling, not considered in [l4j|. is also quite distinct from 
the quenched-in disorder mediated interactions proposed 
in [13| . Since the dipole moment of a single-domain par- 
ticle is proportional to its volume, the effect of interac- 
tion, within a mean-field picture, may be incorporated 
by adding a term proportional to V 2 in the exponent 
of t(V). Thus, even small particles (V2) can now have 
Ti larger than the measurement time. This becomes 
more prominent at lower temperatures. Therefore, the 
blocking temperatures for both particles 1 and 2 are now 
shifted to higher T, thereby causing the wiggles to disap- 
pear. This is consistent with the results of FigQ] which 
show that the memory effects are stronger for the non- 
interacting particles. We conclude then that the unex- 
pected wiggles seen in the cooling and heating cycles of 
M(T) versus T have much less to do with interaction 
effects but more to do with polydispersity of the sample. 
How crucially does the nature of the particle size distri- 
bution function affect the magnetization recovery during 
the zero field heating cycle ? In order to answer this 
question we first quantify the memory effect by defining 
a parameter, 

dM dM 
i? -° ( dT |T = Tj dT' (2) 

where <d(x) is the Heaviside step function. The param- 
eter R measures the positive slope of the M(T) curve 
during zero field heating. We have calculated R using a 
Gaussian size distribution centered at V — Vq and with 
width ,s. Our results for R are shown in Fig[2Id) for a par- 
ticular choice of Vq as a function of s. We observe that R 
increases with the width of the distribution and saturates 
quickly. In this regime, R is almost independent of Vq 
and accordingly, the detailed nature of the distribution. 
We conclude that the memory effects will be best seen in 
samples with a dilute dispersion of particles but a very 
wide (flat) distribution of sizes. Indeed in this limit the 
relaxation is known [l2]| to be prominently dominated by 
magnetic viscosity characterized by a logarithmic relax- 
ation in time. Not surprisingly, a logarithmic relaxation 
has been observed in the experiments of Sun et. oi.[T3| 
although the interpretation offered is different from ours. 
In order to substantiate our interpretation of the M 
vs. T data we have carried out hysteresis measurements 
and thereby coercivity estimation for both the interact- 
ing sample A and non-interacting sample B. The room- 
temperature DC magnetizations versus the applied mag- 
netic field for both samples are shown in FigU3 Clearly, 
for sample B the relaxation times are shorter than the 
measurement time, at 300K. Thus, there is no hystere- 
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FIG. 3: (color on-line) Room temperature M-H curve of sam- 
ples A and B. The hysteretic response of sample A points to 
the presence of strong dipolar interactions. 
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FIG. 4: (color on-line) Coercivity(_f/ c ) as a function of tem- 
perature for the interacting (□) and non-interacting (O) sam- 
ples. The corresponding curves (H c in arbitrary unit) ob- 
tained from our theory assuming a double delta function par- 
ticle size distribution is shown in the inset. 

sis loop and the coercivity (measured by the width along 
the abscissa on the zero-magnetization line) is also zero. 
On the other hand, for sample A, we observe a non-zero 
coercivity even at 300K due to the slowing down of re- 
laxation because of the presence of an additional term 
proportional to V 2 in the exponent of t(V) as mentioned 
above. 

Next we repeat the above measurements down to 4K, us- 
ing a SQUID magnetometer. The coercivity(-ffc) is plot- 
ted as a function of temperature(T), in FigQ] Because 
relaxation slows down for both sample A and sample B, 
H c increases with decrease of T (Fig0J). The coerciv- 
ity of the interacting sample A is larger than that non- 
interacting sample B for temperatures greater than 25K. 
However, at T=25K a surprising crossover is detected, 
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where the coercivity for sample B shoots above that for 
sample A. We suggest that the reason for this behavior 
is that the term H in the exponent of t(V) is replaced by 
H+<5H, where the mean field 5H, arises from interaction: 

5 H =tl v^ { tm±m } (3) 

Kb-1 

where 7 is a parameter that depends on the mean sep- 
aration between the particles. The tanh term augments 
the V 2 term below 25K, making the larger particles relax 
so slowly that they don't respond to H at all. Therefore, 
the larger particles are 'frozen out' from further consider- 
ation, making the mean relaxation time in the interacting 
case even smaller than that for the non-interacting case. 
This somewhat nonintuitive conclusion is further con- 
firmed by our simulated coercivity computation, shown 
in FigQ] (inset). 

To verify our argument further we perform a separate 
set of experiments on both samples A and B as follows. 
We field-cool the samples down to 10K from 300K in the 
presence of H = 100 Oe. At 10K the magnetic field is 
switched off and the relaxation of the magnetization mea- 
sured. We find that the average relaxation time obtained 
by forcing an exponential fit to our data of sample A is 
100 min and that of sample B is 25min. We then heat the 
samples to 300K, and cool it back down to 10K at zero 
magnetic field. At 10K we switch on the magnetic field 
and wait for 2h. The magnetic field is then switched off 
and the magnetization measured. The relaxation time of 
sample B remains 25min. but the relaxation time of sam- 
ple A decreases to 30 min. This result is consistent with 
the reasoning described in the above paragraph. There- 
fore, for the low-temperature interacting system, larger 
particles are rendered magnetically inactive. 
In conclusion, the strong history dependent effects seen in 
magnetization and coercivity measurements in NiFe20<± 
magnetic nanoparticles are interpreted as being due to 
arrested Neel relaxation. Our model is dramatically sim- 
plified by choosing just two volumes of the particles, on 
either side of the 'blocking' limit. Further corrobora- 
tion of the proposed mechanism is achieved by perform- 
ing measurements on an interacting system. Our results 
suggest that either by tuning the interaction (through 
changing inter-particle distance) or by tailoring the par- 
ticle size distribution, these nanosized magnetic systems 
can be put to important application in memory devices. 
In particular, a flat volume distribution can be of great 
utility than a monodispersed distribution with a single 



sharp peak. 
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